Dirac Fermions in Graphite: the State of Art 
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Abstract 

Macroscopic concentration of massless charge carriers with linear conic spectrum - Dirac Fermions (DF) was shown 
in 2004 to exist in Highly Oriented Pyrolytic Graphite (HOPG) and governs its electronic properties. These carriers 
can have the same nature as DF observed in graphite monolayer (Graphene) and let to view HOPG as superposition 
of 2D carbon layers, almost independent electronically. We overview here the recent experimental evidences of 2D 
Dirac Fermions in Graphite and their similarity with carriers in Graphene. 
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1. Do Dirac Fermions present in graphite 

Recent synthesis of graphite monolayer (Graphene) 
[l|2j activated study of several challenging funda- 
mental problems in Solid State Physics related with 
electronic properties of semimetals i.e. systems in 
which several branches of the spectrum cross close 
to Fermi level. It was recognized since a long time 
[3] and confirmed experimentally recently |ri2j that 
peculiarity of electrons in graphite monolayer con- 
sists in the particular linear spectrum of charge 
carriers expanding in vicinity of the corner point K 
of the hexagonal Brillouin zone (BZ), similar to the 
conical dispersion of massless Dirac Fermions (DF) 
in 2D Quantum Electrodynamics. 

Note however that according to the tight-binding 
calculations [4], the linear Dirac spectrum is the 
specific property of Graphene. Already the bilayer 
graphite films have the electronic dispersion that 
is presented as junction of up- and down curved 
parabolas. Although the Dirac-like branches can 
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also appear in the n-layered {n > 2) graphite, their 
relative phase volume is negligible and no manifes- 
tation of DF should be observed. 
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Fig. 1. According to Slonzewski, Weiss and McClure (SWM) 
band model [5], Fermi surface (FS) in graphite is expanded 
along fcz edge K-H of the hexagonal Brillouin zone (BZ) and 
consists of electron and hole pockets. The maximal cross 
section of pockets, seen in quantum-oscillation experiments 
have the massive dispersion of charge carriers whereas Dirac 
Fermions (DF) are concentrated in the vicinity of H-point 
with very small phase volume. 
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At large n the multilayered system recovers 
the properties of bulk graphite and the classical 
Slonzewski, Weiss and McClure (SWM) band model 
[5] should be applicable. According to this model 
the Fermi surface (FS) of graphite is expanded 
around z-directed K-H edge of 3D BZ and have two 
almost compensated electron and hole pockets as 
shown in Fig. [TJ The perpendicular Dirac-like dis- 
persion exists only in the immediate vicinity of the 
point H and therefore the relative concentration of 
DF should be negligibly small. 

Two majority groups of carriers seen by many re- 
searchers in de Haas van Alphen (dHvA) quantum- 
oscillation experiments in 50's-60's [3] (cf. also 
Fig. [2]) were naturally attributed to the maximal 
electron- and hole- cross sections of FS pockets, but 
no definite conclusion about their spectrum (Dirac 
or massive) was done. 

In the beginning of 2004, basing on the phase 
analysis of Shubnikov de Haas (SdH) and dHvA 
oscillations in bulk Highly Oriented Pyrolytic 
Graphite (HOPG) two of us (IL&YK ) discovered 
that one of these groups corresponds to the macro- 
scopic quantity of DF ^ . Evident discrepancy with 
SWM 3D band model concerning the quantity of DF 
in graphite let to the alternative view of graphite as 
the stack of almost independent mono- and multi- 
2D carbon layers. 

Since then, several independent confirmation of 
existence of DF in graphite where done. The objec- 
tive of this communication is to overview the most 
important experimental results that let to judge 
about the nature of charge carriers in graphite. 
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Fig. 2. Fourier spectrum of de Haas van Alphen (dHvA) and 
Shubnikov de Haas quantum oscillations in HOPG reveals 
two types of carriers: normal electrons and Dirac holes. The 
identification of carriers is done by the phase analysis in [6] . 



2. Dirac Fermions from quantum oscillation 

Most of the methods permitting to distinguish be- 
tween normal (massive) carriers and DF are based 
on Landau Level (LL) quantization in magnetic 
field. In the massive case the equidistant LLs 

En ^ {eh/m^c)B{n + 1/2) (1) 

are separated by the gap Eq — ehB /2m±c from 
E = whereas in the Dirac-like case the square root 
dependence takes place: 

En = ±v^2ehBnlc (2) 

and the lowest LL is located exactly at Eq = 0. 
This, in particularly, leads to the difference in Bohr- 
Sommerfeld semiclassical quantization: 

S{Ef) = {n + j)27rh—, (3) 

for which 7 = 1/2 for massive carriers and 7 = in 
the Dirac case [7] . This factor is uniquely related to 
the topological Berry phase = kir (with k = 2j+ 
1 (mod 2)) acquired by a fermion, moving around 
S^Ef) [7j and is manifested experimentally as the 
phase in SdH oscillations of conductivity 

Aa,,(B)~-A(B)cos[2^(:§-7 + <5)] (4) 

or in another quantum oscillations. (FS dependent 
factor \6\ < 1/8) 

Determination of 7 from SdH and dHvA oscilla- 
tions was the base of our phase analysis [5]. Since, 
as was already mentioned, two majority groups of 
carriers are present in graphite (Fig. [2]), we first fil- 
tered the oscillating signal from each of them and 
then measured their phases separately. Moreover, 
comparison of SdH and dHvA experiments let us 
to determine the sign of carriers. Conclusion was 
that, the lower-frequency peak in Fig. [2] corresponds 
to massive electrons whereas the higher-frequency 
peak to Dirac holes. Coincidence of results of SdH 
and dHvA phase analysis shows that existence of 
DF in graphite is the bulk and sample-independent 
effect. 

Contradiction with band theory is not only in the 
predicted by SWM model massive character of both 
groups of carriers [8 (their maximal cross sections 
Se, Sh are far from Dirac point H, see Fig. [1]), but 
also in the inverse order of proportional to Se and 
Sh oscillation frequencies: it was previously assumed 
[3] that Se > Sh- 
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3. Landau Level spectroscopy 

Alternative way to distinguish between different 
types of carriers is the direct determination of the 
energy distribution of quantized LLs. In case of mas- 
sive carriers the inter-LL distance scales linearly 
with B as in Eq. ^ whereas in the Dirac-like case 
as VB as in Eq. 

Low temperature (4K) Scanning Tunnelling 
Spectroscopy (STS) measurements on the surface 
of HOPG in fields up to 12 Tesla were performed 
to test such LLs energy distribution [9]. It was dis- 
covered that the resulting signal reveals the traces 
of both massive and massless LLs, that agrees with 
our phase analysis Unfortunately the natural re- 
striction of STS measurements over the sample 
surface permits also the alternative interpretation 
of the origin of DF as due to the surface effects [5] . 

Another measurements of LLs distribution in 
graphite (both in HOPG and in natural) that are 
free from the mentioned above surface sensitivity 
were done using the Far-Infrared (FIR) magneto- 
transmission spectroscopy |10j in which the ab- 
sorbtion lines correspond to the optical transition 
between different LLs. Two series of LLs with nor- 
mal and Dirac-like quantization in B were found 
and attributed to point K and H of 3D BZ in SWM 
band model. The density of DF was found to be 
smaller then that, extracted from analysis of SdH 
oscillations [5]. 

Although authors of [10] agree that the domi- 
nant transitions in the Dirac series of absorbtion 
lines have their counterpart in the FIR spectra of 
graphene [TT] , they believe that the 3D band model 
is more adequate to describe the situation because 
of observation of the weak series of additional transi- 
tions, forbidden in pure 2D case. Note however that 
in this case one more series of FIR transitions, cor- 
responded to the maximal cross section of the hole 
FS pocket Sh (Fig.[T]) with massive LL quantization 
should be observed. 

4. Angle Resolved Photoemission 

Coexisting of massless DF with quasiparticles 
with finite effective mass in graphite was also proven 
by Angle Resolved Photoemission Spectroscopy 
(ARPES) [H]. This method is the most direct and 
unique tool to measure the electronic structure with 
both energy and momentum information. The elec- 
tronic spectrum was scanned in the inverse space 



of BZ in vicinity of points K and H in the k _L 
direction. The electronic dispersion was shown does 
correspond to the massive carriers in K-point and 
to DF in H-point as was predicted by SWM band 
model. 

Although ARPES measurements convincingly 
demonstrate the presence of DF in graphite, no fur- 
ther information about their quantity and macro- 
scopic properties can be extracted. To know their 
concentration and to check the validity of the as- 
sumed in 12] 3D band model the further scans in 
k _L k^ direction along whole K-H edge are required. 

5. 2D Transport and Quantum Hall Effect 

To give more insight in favor of 2D independent 
layered structure of graphite we refer on another 
group of experiments. Note first that because of 
the extremely high out- to in- plane transport 
anisotropy in HOPG: pout/pin ^ 5-10'' 13J (instead 
of 300 as was measured for some graphite crystals 
[3] ) it is unlikely that out-of-plane transport can be 
explained within the inter-plane hopping of SWM 
model. Moreover, the characteristic property of 
2D conducting systems - the Quantum Hall Effect 
(QHE) has been discovered in 2003 in HOPG [13]. 

Careful analysis of QHE staircase in HOPG [H] 
shows that it contains the precursors of both QHE 
discovered in graphene |ll2j and in bilayer graphite 
[15] , having the semi- integer (Dirac) |16|17j and in- 
teger (massive) 1181 LL quantization correspond- 
ingly. Order and appearance of normal and Dirac 
steps in HOPG depends on the sample [19 and 
sometimes coincides with QHE staircase observed in 
few-layer graphite sample 20J. Such behavior con- 
firm the hypothesis that bulk graphite is a system 
randomly composed of single- and few- carbon lay- 
ers. 

6. Raman fingerprint of Graphene in 
graphite 

Next, we present the results of our recent micro- 
Raman testing of the best (in 2D transport sense) 
HOPG samples. (The details will be published 
elsewhere.) Raman study of the second order peak 
"2D", which is provided by specific to Dirac spec- 
trum double-resonance electron transitions is 
known to be the best tool to distinguish the mono- 
and multilayered graphite films [22] . 
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Fig. 3. Double resonance 2D micro-Raman band of HOPG 
graphite (two blue plots in the middle) demonstrates the 
superposition of the features typical for bulk natural graphite 
(upper plot with left-shoulder profile) and for the monolayer 
graphene (lower single-peak plot) |22) . The relative weight 
of these features depends on the position of the laser spot 
on the surface of HOPG. The typical Raman spectrum of 
graphite is shown in the upper-left corner 

Shown in Fig. [3] the upper left-shoulder-peak 
"2D" in bulk graphite is indeed very different 
from lower, specific only to monolayer the single- 
Lorentzian peak of graphene [22 . Two middle plots 
correspond to our micro-Raman testing of HOPG 
with different position of few-micron laser spot on 
the surface. 

Surprisingly we found that depending on the spot 
position, the classical graphite peak can demon- 
strate the growing fingerprint feature of graphene. 
Because of several microns laser penetration length, 
such effect can be explain only by existence of 
macroscopic number of independent graphene lay- 
ers in bulk of HOPG. Since this conclusion is coher- 
ent with other considered above experimental facts 
we can express even more general hypothesis that 
left shoulder in peak "2D" in graphite is provided 
in total by the presence of graphene monolayers in 
the bulk. 

7. Conclusion 

Retrospective overview of the experimental data 
let to believe that macroscopic concentration of DF 
does exist in the bulk of graphite. Since this con- 
clusion contradicts to classical SWM band model, 
further experimental and theoretical efforts should 
be done to clarify their nature. It is of primary im- 
portance to understand, whether the modified band 
theory should be constructed to explain the exis- 



tence of DF or indeed, graphite can be viewed as 
the system of independent mono- and multi carbon 
layers as proposed in the present article. 
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